The importance of being discrete - life always wins on the surface 
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Many systems in chemistry, biology, finance and social sciences present emerging features which 
are not easy to guess from the elementary interactions of their microscopic individual components. 

In the past, the macroscopic behavior of such systems was modeled by assuming that the collective 
dynamics of microscopic components can be effectively described collectively by equations acting on 
spatially continuous density distributions. It turns out that quite contrary, taking into account the 
actual individual/discrete character of the microscopic components of these systems is crucial for 
explaining their macroscopic behavior. 

In fact, we find that in conditions in which the continuum approach would predict the extinction 
of all the population (respectively the vanishing of the invested capital or of the concentration of a 
chemical substance, etc), the microscopic granularity insures the emergence of macroscopic localized 
sub-populations with collective adaptive properties which allow their survival and development. In 
particular it is found that in 2 dimensions "life" (the localized proliferating phase) always prevails. 



In addition to physics, an increasing range of sciences: 
chemistry, biology, ecology, finance, urban and social 
planning have passed in the last century to quantitative 
mathematical methods. 

Along with the obvious benefits, it turns out that the 
traditional differential equations approach has brought 
some fallacy into their study. 

We present here a very simple generic model which con- 
tains proliferating (and dying) individuals and we show 
that in reality it behaves very differently than its rep- 
resentation in terms of continuum density distributions: 
In conditions in which the continuum equations predict 
the population extinction, the individuals self-organize 
in spatio-temporally localized adaptive patches which in- 
sure their survival and development. 

This phenomenon admits multiple interpretations in 
various fields: 

-if the individuals are interpreted as interacting 
molecules, the resulting chemical system emerges spatial 
patches of high density which evolve adaptively in a way 
similar with the first self-sustaing systems which might 
have anticipated living cells. 

- if the individuals are the carriers of specific genotypes 
represented in the genetic space, the patches can be iden- 
tified with species, which rather than becoming extinct, 
evolve between various genomes (locations in the genetic 
space) by abandoning regions of low viability in favor of 
more viable regions. This adaptive speciation behavior 
emerges in spite of the total randomness we assume for 
the individuals motions in the genetic space (mutations). 

- interpreted as financial traders, the individuals de- 
velop a "herding" behavior in spite of the fact that we 
do not introduce communication or interaction between 
them. This leads to the flourishing of markets which the 
continuum analysis would doom to extinction. 

All these phenomena have in common the emergence 
of large, macroscopic structures from apparently uni- 
form backgroundO due to the amplification of small, mi- 
croscopic fluctuations which originate in the individual- 



ized character of the elementary components of the sys- 
tem. This mechanism insures in particular that on large 
enough 2 dimensional surfaces, even if the average growth 
rate is negative (due to very large death rate), adaptive 
structures always emerge and flourish. 

Imagine an area inhabited by a population of eternal 
agents A which are spread out uniformly with average 
density tla and move around randomly, with diffusion 
coefficient Da- Imagine now a race of mortals, B, which 
are also spread over this area, with initial uniform den- 
sity n#(0). The B agents die at a constant rate, /i, 
(B — — > 0) and proliferate (divide) when they meet the 

"catalysator" , A, with rate A (B + A B + B + A). 
The B's are diffusive, hopping at the rate D B . What will 
happen? 

The naive lore based on macroscopic continuity as- 
sumptions will predict that A reaches a spatially homo- 
geneous distribution, ua(x) = ua, while the B time vari- 
ation ^f- is represented by the linear differential partial 
differential equation: 



dn B 
dt 



D B V 2 n B + (Xn A - n)n B 



(1) 



The first term represents the uniformization effect of B 
diffusion while the iiB term represents the fact that a cer- 
tain fraction of B's die per unit time. The crucial term 
XnATi B represents the proliferation of B's in the presence 
of the "life giving" A's. Note that the equation is linear 
in n B and that for initial spatially uniform tla and n B 
distributions it has the time exponential solution 

n B (t)=n B (0)e^ XnA -^. (2) 

In particular Eq. (||), predicts that if the macroscopic 
proliferation rate Xua is lower than the death rate /i, the 
B population will uniformly decrease to extinction. 

Using microscopic representation techniquesH, one 
finds that populations of discrete proliferating agents are 
much more resilient than one would first guess based on 
macroscopic or continuum (PDE) treatment (Fig. 1). 
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FIG. 1. The lower half of the figure shows the exponential 
growth of the average B population as a function of time, in 
the actual simulation (solid blue line) compared to the expo- 
nential decrease predicted by the continuum approximation 
(dashed red line). < A > is the average number of A reac- 
tants per site. Other symbols are defined in the text. The 
snapshots in the upper half show the spatial configuration 
of A and B react ans. B react ans are seen to be localized 
in islands (notice that what is plotted is the logarithm of B 
concentrations, thus localization is stronger than would first 
appear) . 

The study of diffusion limited reactions! have already 
shown in the past deviations from the continuum the- 
ory due to the quantized nature of the react ants. In the 
present case the effect is even more dramatic: it consti- 
tutes the difference between life and death: the contin- 
uum approach predicts extinction while the direct sim- 
ulation uncovers the emergence of a thriving, adaptive, 
developing system (Fig 1). 

In order to understand the " source of life" in this sys- 
tem one has to concentrate on the microscopic conditions 
around the individual A agents, rather then looking at 
the local average growth rate Xua — 

Fig. 2 represents the evolution of the B cloud follow- 
ing a single A agent as it jumps around randomly. The 
B concentration is shown to trace the A as it performs 
a random walk. Clearly, the colony does not decay to 
extinction ; instead, it seems to trace the A, "trying" to 
keep its center of mass at its location. As shown in the 
inset, each jump of the A is followed by a momentary de- 
crease in the height of the B concentration, however, due 
to the multiplicative process there is an overall increase. 




FIG. 2. The profile of a B island as a function of time as it 
follows the random motion of an A agent. The cross-section 
of the island is taken through the current location of the A 
agent. The inset shows the time evolution of the height of the 
B concentration at the point at which A is currently located 
(solid blue line). The B colony is seen to grow, although the 
average growth rate over the entire space is negative {ua is 
extrenly low since there is only one in the whole simulation 
space, thus Xua — /jl ~ — /x). The dashed red line shows the 
and exponential growth with coefficient eo — 2(1Dak where 
k, — Log e {\/ Db), is the slope of the island (this slope is ex- 
hibted in the main graph, and can be derived from a simple 
approximate calculation), eo is A — \i — 2(1Db derived similarly. 

Let us consider first the simplest situation of a single A 
agent jumping randomly (with a diffusion coefficent Da) 
between the locations in an infinite d-dimensional space. 
In between A lumps the B density at the A location grows 
exponentially! as n B (t) ~ n B (^)e^-^- 2dDB ^ . Where A, 
fi and 2dDs stand for the proliferation, death, and the 
loss due to diffusion, respectively. The estimation is made 
by neglecting the flow of B's from neighbouring site to 
the A site, this is justified when the B concentration in 
the neighbouring sites is much lower than on the A site. 
In the same limit, the ratio between the height of the B 
density at the A location and the height of the B density 
on a neighboring site is easily estimated: X/Db- Conse- 
quently, each A jump corresponds to a sudden downwards 
jump by a factor of X/Db in the height of the B hill . As 
there are in average 2(1Da such jumps per unit time, the 
net effect of proliferation, diffusion and death, gives the 
B concentration at the A site as a function of time: 

n B (t) = nB (o)e^ x -^- 2dDB - 2dDALo9eWDB))t (3) 

The approximation is in good agreement with the sim- 
ulation shown in Fig. 2. The slope of the island, on 
a log scale, is indeed seen to be Log e (Xj Db\ the the 
time dependence of the height of the B island in between 
A jumps is indeed given approximately by an exponent, 
(X — \i — 2dD b)^. Consequently the dashed red line (in the 



x The above analysis turns void if Db — 0, where the spa- 
tial dimensions of the island do not grow at all. We do not 
consider this singular case in this paper. 
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inset) which represents Eq. 3 follows closely the actual 
growth seen in the simulation (blue line). The difference 
between the theory and simulation is mainly due to cases 
where two or more A jumps follow each other rapidly, in 
this case the island's shape does not stabilize before an- 
other A jump is made, these rather rare events, modify 
somewhat the actualy result 1 . 

One may ask what is the situation in the case, when 
single colonies are unstable (i.e where the exponent in 
Eq. 3 is negative). One possibility is that in such a sit- 
uation the continuum approximation is valid and the B 
concentration decays to zero. Another possibility is that, 
although single isolated colonies are unstable, global ef- 
fects such as islands growing, joining and splitting give 
us back the survival feature. In particular, since large 
colonies are more stable than small colonies, one may ex- 
pect the typical size of an "active" colony to grow with 
time. This behavior is demonstrated in Fig. 3 which 
show the active clusters in a two dimensional system de- 
veloping in time. Evidently, the small clusters either de- 
cay or merge into larger and larger clusters. 





FIG. 3. The spatial distribution of B's for different times 
(earlier times to the left), where concentrations greater than 
10 are colored red and concentrations below 10 are color 
coded. 

The arguments and computer experiments mentioned 
above show convincingly that the individuals based life 
is much more resilient than a hypothetical life density 
spread continuously across spatial regions. 

These conclusions should suffice to induce profession- 
als in biology, finance and social sciences to consider seri- 
ously the regime in which their systems are naively non- 
viable (decay to extinction) when viewed macroscopically 
but perfectly viable in reality (and wheii, simulated cor- 
rectly at the microscopic individual levelu). 

In particular, most of the species in nature could be in 
this regime: negative naive average macroscopic growth 
rate but actual survival and proliferation. Similarly, mar- 
kets which might look unappealing when averaging over 
the various investing possibilities might prove lucrative 
enough (at least for the lucky investors which hit profit 
opportunities A) as to maintain them in the competi- 
tive range. In fact this line of thought might provide 
an explanation to the emergence of life from the ran- 
dom chemistry of its component molecules in spite of the 



formal extreme improbability of the event. Equally it 
might explain the paradoxes in finance between the effi- 
cient market hypothesis (absence of systematic profit op- 
portunities in equilibrium markets) and the actual profits 
which investors extract daily from the market. 

In order to obtain a more rigorous bound on the param- 
eter range in which life overcomes the gloom prognosis 
of the macroscopic analysis we used the renormalization 
group (RG) analysis which indicates that on large enough 
surfaces, life always wins. For higher dimensions, the do- 
minion of life still extends to arbitrary low ua densities, 
but a minimal finite A value is required. 

In RG, the collective behavior of the system is identi- 
fied by integrating out the small length scale, short time 
fluctuations, leaving us with an effective theory for the 
large scale objects. Here, these are the large, stable is- 
lands shown in Fig. 3. The new, effective theory is char- 
acterized by renormalized coupling constants, i.e., modi- 
fied numerical values of the effective rates (growth rate, 
dearth rate, hopping etc.) on large length scale. The 
process of decimating small fluctuations is then iterated 
again and again, giving us flow line which reflect the evo- 
lution of the effective values of the coupling constants as 
one integrates larger and larger scales I. 
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Fig. 4 shows the flow lines of m (m = \i — Xua) and A 
due to the iteration of the decimation process 3 . The flow 
is given by the equations: 



dm 
~dF 



2m 



dX 
~dl 



A[2-d + 



A 2 n A 
2nD 

— ] 



(4) 



For d < 2 we see that for large length and time scales 
(that is, after many iterations of the decimation process), 
A grows without limit while m eventally becomes nega- 
tive. This implies that on the large scale, the system 
actually behaves as if Xua > and life always wins. 



2 The details of this RG analysis, which involves the pre- 
sentation of the exact Master equation of the process as a 
field integral and the e-expansion around the critical dimen- 
sion d c — 2 are out of the scope of this report and will be 
presented elsewhere. 

3 D = Da + Db is the effective diffusion constant 
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FIG. 4. Lower panel shows flow lines for d > 2. Shaded 
region flows to negative mass ("life"). Upper panel shows flow 
lines for d < 2, the whole parameter space flows to negative 
mass. 

In higher dimensions (d > 2) Fig. 4a indicates a dy- 
namical phase transition where for part of the parameter 
space the system flows to negative m (life) and for an- 
other part the system flows to positive m (death). 

It should be noted that the flow portrayed in Fig. 4 
is associated with larger and larger length scales: For a 
finite system, the flows should be truncated and the size 



of the system may be crucial: simulations with parame- 
ters identical to that of Fig. 3, lead to extinction when 
carried out on a system size 4 times smaller. 

In conclusion, our results suggest that the dimension- 
ality of the system and its size are crucial features for 
its capability to emerge and sustain life. This may ex- 
plain the fact that most of the ecological systems are 
two dimensional. Reinterpreting in the genome space, 
the present results provide the conceptual link between 
the atomized structure of the life building blocks and the 
explosive Darwinian tandem, noise + proliferation. 
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